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Abstract 

We present the results of a next-to-leading order calculation of Wt production, including the 
decays of both the top quark and the W boson. The effects of radiation in the decay of the top 
quark are also included. The separation of diagrams which appear in the real corrections, into 
singly- and doubly-resonant contributions, is performed using a 6-jet veto which is motivated by 
the use of the bottom quark distribution function. We find that, for a choice of scale which is 
suitable for this approach, the QCD corrections are very mild and only change the cross section by 
up to 10% at the LHC, depending on the severity of the 6-jet veto. When further cuts are applied, 
applicable for a Higgs boson search in the H — ► WW* channel, we find that the radiative effects 
greatly decrease the number of background events expected from this process. In addition, the 
shapes of relevant distributions can be significantly changed at next-to-leading order. 
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I. INTRODUCTION 



At the LHC, the top quark will be produced copiously in many channels. As well as the 
it pair production channel, the top quark may be produced singly in association with other 
particles. The rates for these processes will be sufficient to both study the properties of 
the top quark in detail and to provide a significant source of background events for other 
analyses 

In this paper we will discuss the calculation of the next-to-leading order (NLO) corrections 
to the production of a single top quark in association with a W boson. This calculation 
has been included in the general purpose NLO program MCFM . The lowest order 

process which we consider is, 

b + g -> W- + t 

Ui/ + e + + 6 (1) 



v 



so that the leptonic decays of the W~ and of the top quark are included. We note that both 
at the Tevatron and the LHC, the rate for the charge-conjugate process involving a W + 
and a t quark is identical j^. This process has previously been considered extensively at 
leading order 0,0,3]. However, it is only at next-to- leading order that we obtain accurate 
predictions of event rates which are sensitive to the structure of jets in the final state. Such 
NLO calculations have so far been available only for the case where the decays of the heavy 
quark and W boson are not included 0, j|. 

We have extended these predictions to include not only the full spin correlations in the 
decays of the W boson and the top quark, but also to include the effects of gluon radiation 
in the top quark decay. This is achieved using the same method that has previously been 



and which is briefly described, together 



applied to other single top production channels 
with other details of the calculation, in Section HTl 

At next-to-leading order some of the contributions representing the emission of an addi- 
tional parton require special attention. One finds that the corrections involving two gluons 
in the initial state contain diagrams that would normally be assigned to the lowest order 
calculation of the doubly resonant ti production process. We discuss our treatment of this 
complication in Section ITTTl 

The results of our calculation are presented in SectionllVl We discuss the NLO corrections 



at the Tevatron and the LHC, comparing our findings with those obtained previously in the 
literature. We also provide updated predictions obtained using the latest experimental 
inputs and examine the effect of including gluon radiation in the decay of the top quark. 

Section |3 contains a study of our results in the context of the search for an intermediate 
mass Higgs boson at the LHC. In this channel the Higgs decays via WW*, with the final 
state containing leptons and missing transverse momentum. Since the Higgs mass cannot 
then be reconstructed, theoretical input as accurate as possible is imperative. To that end, in 
this section we apply realistic acceptance and search cuts to all the final state particles, then 
compare the effect of the NLO corrections with the more inclusive results already presented. 



II. CALCULATIONAL DETAILS 



To evaluate the matrix elements for the production and decay of a H^-to p sy stem, we 
follow the same strategy as in a previous calculation of single top production 10( , which is 
based on two approximations. The first is that the top quark is produced and decays exactly 
on-shell, motivated by the fact that diagrams without an on shell top quark are suppressed 
by a factor of T t /m t Pd 1%. This enables a division of the process into production and decay 
stages, with the further approximation that the interference between radiation in the two 
stages is neglected. On general grounds the contribution of these interference terms can be 
shown to be of the order of ot s T t /m t due to the large difference between the characteristic 
time scales of the production (m^ 1 ) and decay (IT 1 ) stages. More technical details and 
further references can be found in Ref. jlfl| . 

The tree level amplitude is represented in Fig. Q where the on-shell top propagator is 
denoted by two short lines. Due to the weak vertices the only two non-vanishing helicity 
amplitudes correspond to the two polarizations of the gluon, with all massless fermions 
left-handed. In our expressions for the amplitudes we take all momenta outgoing, restoring 
incoming momenta in our implementation by performing the proper analytical continuation 
of the massless spinor products. We write the amplitudes in terms of the momenta of the 
decay products, with labels as follows, 

bi + g -> W- + t 

— > Vf + Bf + bf (2) 

-> ej + Vi 
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FIG. 1: The two tree level diagrams for the Wt process. 



and also use t to represent the momentum of the top quark, so that t = i/f + ef + bf. 
Both bottom quarks are treated as massless particles in our approach. The two tree-level 
amplitudes are then given by, 

A- = t—tt (m 2 t (gei) [b i ef]-(g\t\e f )(e i u i ) [&ii/ f ]), 
L 9 Oi \ 

f ( ' HI _ \ / _ H I _ \ _ 2 / L _ \ r _ „ 1 , O _ , I _ \j. I _ \ [P ^ 



^4+ = (< &* l^le/Xeilt | (/) -m t (feiei) [^e/] + 2g • t ( e< 1 1 j e f ) ^ J , (3) 

where the overall factor / is, 

/ = g s gjT a {u f b f ) [bjVj] _ 

y2 g • t Wprop Wprgp t prop 

In these formulae P prop = P 2 — Mp + iMpTp, T a represents one of the eight SU(3) color 
generators and the spinor products with massless four-momenta pi and pj are defined as 
usual: 

(PiPj) = (Pi ~ \ Pj+), (5) 
[PiPj] = (Pi + \Pj~), (6) 
(Pi\Pk\Pj) = (Pi~ \fa\Pj-}- (7) 

As explained above, the virtual corrections to this process can be divided into production 
stage corrections (represented in Fig. |2J) and decay stage ones (shown in Fig. The same 
division applies to the real corrections, which are depicted in Figs. |U and El respectively. By 
producing the top quark strictly on shell we are assured that the diagrams in Figs. 0] and El 
are separately gauge invariant. 
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FIG. 2: One loop diagrams representing virtual corrections in the production stage. 



The cancellation of the soft and collinear divergences between real and virtual contri- 

n 

butions has been implemented through the subtraction method In particular, for the 
contribution from real radiation in the production stage we have adopted an extension of 



the dipo 

state 



e subtraction scheme 



12. 



1^ which handles the case of massive quarks in the final 
14j . We have used a generalization of this method, where one can use a tuneable 
parameter in order to have better control over the size of the subtraction. Further details 
and formulae are contained in Ref. and Appendix [X] For the case of real radiation in 
the decay we have used a further extension of this method, as in Ref. [l^, which ensures 



FIG. 3: One loop diagrams contributing to the calculation of virtual radiation in the decay stage. 




FIG. 4: Diagrams representing the emission of real radiation in the production stage. In the 
calculation, the additional gluon must also be crossed into the initial state and diagrams containing 
two quark lines (not shown) are also included. 

that the top quark and the W boson remain on-shell when the subtraction is performed. 

The virtual corrections in the decay stage amount to the study of the decay of an on-shell 
top quark. These amplitudes have been known for a long time [JJj and we do not report 
them here. To evaluate the contribution of the virtual radiation in the production stage 
(Fig. |2I) we start from the amplitudes where the top quark is produced on-shell without 
decaying and calculate amplitudes for the two polarization states of the top quark. This 
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FIG. 5: Diagrams contributing to real radiation in the decay of the top quark. 



is achieved by writing the spinors in terms of an auxiliary massless four-vector to in the 
following way: 



The vector to is constructed by forming a linear combination of t and g, 



The full result, where the decay of the top quark is included, can then be obtained by 
combining these amplitudes with the ones for the decay t — > Wb, calculated in the same 
way. This is possible since the intermediate top quark propagator is recovered via the 
identity, 



Performing the calculation in this factorized way has a number of advantages. First, useful 
consistency checks can be performed using the amplitudes without the top quark decay. 
Second, by replacing the top quark mass appropriately, these amplitudes can be used to 
study other processes where this decay is not relevant. 

The diagrams of Fig. El have been written and algebraically manipulated using FORM, 
after we have used the background field gauge to contain the number of terms generated 






(10) 
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by the three gluon vertex. We deal with infrared and ultraviolet singularities by using 
dimensional regularization in the four dimensional helicity scheme and use the method of 



Ref. 17] to write the amplitudes in terms of traces. The appearance of 75 in the weak vertex 
is then handled by cyclically rotating these traces so that 75 appears at the beginning of each 
trace, before performing the contraction of Lorentz indices. Using this prescription we have 
checked that the Ward identity for the weak current is satisfied exactly and no additional 
counter terms are required. 

Finally, we are left with box vector integrals and triangle rank 2 tensor integrals. Us- 
ing Passarino-Veltman n-dimensional decomposition, we obtain a result in terms of scalar 
integrals. Due to the nature of our approach, other spurious divergences are still present 
at this stage. Individual terms in the result appear to be divergent as factors in the de- 
nominator approach zero. However, in this limit, a combination of such apparently-singular 
terms is finite. By collecting all such terms over a common denominator, one can identify 
new functions that are well-behaved in these limits. These are combinations of rational and 
logarithmic functions, as in Ref. Following this procedure, we are able to refine our 

first result considerably. To show the level of simplification that we have reached, one of the 
amplitudes is reported in Appendix El The others are slightly larger in size and we do not 
reproduce them here. They will be available, together with the rest of our calculation, as 
part of the next release of the Monte Carlo program MCFM. 



III. SEPARATION OF Wt AND tt DIAGRAMS 



When calculating the real radiation corrections, all appropriate crossings of the diagrams 
shown in Fig. 0] should be included. Some of the crossings, in which the additional parton 
is a b quark in the final state, are particularly problematic. These diagrams are shown 
in Fig. IH1 and involve gluon-gluon and same flavour quark-antiquark initial states, which 
are important at the LHC and Tevatron respectively. All these diagrams produce a final 
state consisting of a W~, an on-shell top quark and a b quark. However, the diagrams in 
panel (b) contain a resonant i propagator and represent the production of a ti pair with the 
subsequent decay of the i into the W~ and b quark. As such, the contribution from these 
diagrams when integrated over the total available phase space can be much larger than the 
lowest order Wt cross section (an order of magnitude at the LHC). In order to disentangle 



S 
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(a) (6) 

FIG. 6: Real corrections to W~t production which involve an additional b quark. The double 
bars indicate the on-shell top quark which subsequently decays into W + b. Diagrams obtained by 
interchanging two gluons are not shown. The 3 representative diagrams in the right panel (b) 
contain a resonant i propagator, while those on the left (a) do not. 



these two processes, two methods have been outlined in the literature. 

The first involves making a cut on the invariant mass of the W~b system to prevent 
the t propagator from becoming resonant p. The second method instead subtracts the 



contribution from the resonant diagrams so that no on-shell piece remains [6] . A comparison 
of these two approaches 2] shows that the methods yield the same total cross section when 
a mass window of 15I\ w 25 GeV either side of the top mass is chosen. However, these 
methods do not lend themselves to a Monte Carlo implementation where one wishes to 
study distributions of final state particles as well as total cross sections. Therefore we shall 
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adopt neither of these prescriptions but instead follow a procedure motivated by our use of 
the bottom quark PDF. 

In the 6-PDF approach, the b quark distribution function is derived perturbatively from 
a collinear g —>■ bb splitting that occurs in the initial state. It implicitly includes all split- 
tings up to a pt of the 6-quark equal to the factorization scale, fip. This means that the 
contribution from the corresponding gg — > Wtb diagrams (contained in panel (a) of Fig. |)J) 
has already been included in the lowest order calculation. Therefore the net contribution 
from these diagrams, including appropriate counter-terms and integrating over all b quark 
transverse momenta up to fiF, should be approximately zero. For a suitable choice of \if we 
have checked that this is indeed the case. The choice of fip is made such that the collinear 
approximation used in deriving the 6-PDF is accurate, which for this process implies that 

n 

< {m w + m t )/A « 65 GeV |39|. 

There is also a contribution from the diagrams in panel (b) of Fig. 13 when the p? of the 
b quark is small, p T < //p. However, these diagrams simply rep resent the doubly- resonant 
tt production process which is currently known up to NLO [40]. Therefore it is preferable 
to separate this contribution from the 'genuine' NLO corrections to the Wt process. As 
we shall show later, although the contribution from the tt diagrams in this region of phase 
space is rather small in relation to the total cross section, it is still competitive with the Wt 
result. This suppression means that the interference effects between the two sets of diagrams 
is very small when using a 6-jet veto, in contrast to the case when p\ is unconstrained Q|. 

When a b quark is observed with a px above /ip then our description of the final state is 
a lowest order one. The contribution from the doubly-resonant diagrams dominates and, as 
above, a better prediction would be obtained by using the tt process. Alternatively, one could 
use a calculation including all the diagrams for gg — > tWb, including the t — > bW — > biv 
decay (retaining the b quark mass) and also finite width effects (2]]]. However, currently this 
study would be limited to leading order in a s only. 

To summarize, we shall perform our calculation of the Wt process by applying a veto 
on the pt of the additional b quark that appears at next-to-leading order. This aids the 
separation of this process from doubly-resonant ti production. When applying this veto, 
one should choose the factorization scale equal to (or at least of the same order as) the 
maximum p? of the b quark that is allowed, p T veto . This choice respects the approximations 
that were originally used to define the b quark PDF. For fip 7^ p T veto and for less inclusive 
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quantities, the contribution from the gg — > W^tfr diagrams is calculated by simply omitting 
the doubly-resonant diagrams [41]. The result for this piece remains at the level of a few 
percent of the lowest order cross section. 



IV. RESULTS 

Before discussing the effect of including radiation in the decay of the top quark, we will 
first consider just the Wt total cross section in order to discuss some features of our approach 
and to compare our results with those available in the literature. 

A. Comparison with no top quark decay 

The NLO corrections to the total Wt cross section, where no decay of the top quark is 
included, were previously presented in Ref. jjj. For the sake of comparison, in this section 
we will adopt the parameters therein as closely as possible. In particular, we choose m t = 
175 GeV and the CTEQ5 set of parton distribution functions |42j. The other electroweak 
parameters that enter our calculation are chosen to be, 

M w = 80.419 GeV, g w = 0.4267. (11) 

We perform our comparison at the LHC and examine the dependence of our results on 
the common renormalization and factorization scale /x. As we have already argued, the 
6-PDF approach is most well-motivated when choosing a value of /x less than about 65 GeV. 
However, such a value is much smaller than the more typical choice, /x = m t + mw which is 
the central value chosen in Ref. Q. Therefore, for the sake of illustration, we choose to study 
the scale dependence over a large range from /x = 25 GeV up to /x = m t + mw = 255 GeV. 
As we have discussed above, we limit the px of the b quark that appears at next-to-leading 
order to have a maximum value p T veto , which we choose here to be 50 GeV. 

Our results at LO and NLO are shown in Fig. d We see that the dependence of the lowest 
order curve on a common scale choice is already remarkably small. At next-to-leading order 
we see that this is improved still further, with the cross section varying by about 3% over 
the range of scales shown in the figure. 

Comparing with Zhu jjj, we find a number of differences. Even at lowest order our result 
for /x = m t + mw is higher and furthermore, the dependence of the result on this scale 
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FIG. 7: Scale dependence of the cross sections for W t production at the LHC for rat = 175 GeV. 
The branching ratios for the decay of the top quark and the W's are not included. Cross sections 
are evaluated with CTEQ5L1 (a s (M z ) = 0.127) and CTEQ5M1 (a s (M z ) = 0.118) PDFs 22]. We 
choose a single common renormalization and factorization scale, fi. The lowest order cross section 
is the dashed curve, whilst the NLO one - calculated with 

p&veto = 50 GeV -i s solid. 
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appears less mild (c.f. Fig. 2 of Ref. |9(). However, we have checked the lowest order results 

] and found good agreement. 



of our program against those obtained with MadEvent 
At next-to-leading order we also find a slightly different result, lower and with a stronger 
dependence on the scale. In this case, we expect some discrepancy due to our different 
method of handling the gg — > Wtb contribution. We note that the combination of the b- 
jet veto and our preferred choice of a much smaller scale, leads to a next-to-leading order 
cross section that is about 15% smaller than that found by Zhu. In addition, the i^-factor, 
defined as the ratio of NLO and LO cross sections, is much smaller and in the range 1.2-1.3 
depending on the scale choice. 

Finally, we consider the dependence of our result on the choice for p£ veto . In Fig. El 
we show the scale dependence for three different choices of the veto threshold. We have 
concentrated on the region of smaller scales and varied the common scales by a factor of two 
about the central value /j,q = p|. veto . One sees that, within this window, the scale dependence 
of the next-to-leading order calculation is again somewhat smaller than that found at LO 
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FIG. 8: Scale dependence of the cross sections for W~t production at the LHC, for three different 
choices of veto = /xo- From bottom to top, the solid curves represent values of 30, 40 and 50 GeV. 
The lowest order dashed curve (calculated with /io = 50 GeV) is also shown for reference. 



and improves as the veto threshold is raised. Compared to a threshold of 50 GeV, the cross 
section decreases by about 15% when lowering the veto to 30 GeV. This substantially redues 
the effect of the next-to-leading order corrections on the cross section, leaving a i^-factor 
close to unity for our central scale choice. 



B. Updated results and radiation in the decay 

First, we repeat the calculation of the total Wt cross section, but using the most recent 



determination of the top quark mass 24], which yields m t = 178 GeV. We also use recent 
PDF sets from the MRST and CTEQ groups. 

Our predictions for the Tevatron and the LHC are shown in Table EJ where we have used 
pb veto _ ejQ Q e y anc [ factorization and renormalization scales also equal to this value. It is 
clear from the quoted cross sections that this process is of little phenomenological relevance 
at the Fermilab collider, although we include the result here for completeness. The NLO 
corrections at the Tevatron increase the cross section by a factor of 1.35 when using the 
CTEQ PDF set, but decrease it by a little under 10% for the MRST parametrization. The 
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TABLE I: LO and NLO cross sections for single top-quark production in association with a W~ 
at the Tevatron and LHC, for rat = 178 GeV. The branching ratios for the decays of the top quark 
and the W~ are not included. Cross sections are evaluated with CTEQ6L1 (a s (Mz) = 0.130) and 
CTEQ6M Ja s (M z ) = 0.118) PDFs 3 and also with the MRST2002 NLO PDF set (a s {M z ) = 
0.1197) [2a]. The renormalization and factorization scales are set to the b jet veto threshold of 
50 GeV. The errors represent Monte Carlo statistics only. 



Collider, yfs 


PDF 


<?lo [pb] 


vnlo [pb] 


pp, 1.96 TeV 


CTEQ6 


0.04796 


0.06458 ± 0.0001 


pp, 1.96 TeV 


MRST2002 


0.08083 


0.07414 ± 0.0001 


pp, 14 TeV 


CTEQ6 


29.41 


32.10 ± 0.03 


pp, 14 TeV 


MRST2002 


31.08 


34.49 ± 0.03 



NLO results still differ by about 15%, reflecting the considerable uncertainty in the gluon 
distribution at large x. At the LHC the effect of the NLO corrections is much smaller, 
increasing the cross section by about 10% in both cases, and the predictions from the two 
PDF sets show much better agreement. 

We now turn to a fuller description of the final state, where all the leptonic decays of 
the top quark and the W bosons are included, as in Eq. (JTJ). The decays are included using 
Breit-Wigner propagators with widths, 

T w = 2.06 GeV, = 1.651 GeV, (12) 

and we now investigate the effects of the inclusion of gluon radiation in the decay of the top 
quark. Although the inclusion of this radiation should not change the total cross section, 
a difference is expected when working at a fixed order of perturbation theory. Our results 
are summarized in Table |H] We have used the CTEQ6M PDF set for all the cross sections 
in this table, so that the effect of including radiation in the decay can be understood more 
easily. When including radiation in the decay, the total cross section should change by an 
amount that is formally of higher order in % and is given by [loj | . 

( r t L ° \ 

O Bt^bve+X — CrBt^bue = ( -^NLO ~ M ( a Bt^bue ~ &()Bt^bve) ■ (13) 
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TABLE II: Comparison of LO and NLO cross sections for W~t production at the Tevatron and 
LHC, with leptonic decays of both the W~ and the top quark. The NLO calculation is per- 
formed both without including QCD effects in the decay (aBt^bue) and also when it is included 
(aBt^bue+x)- The top quark mass is mt = 178 GeV and cross sections are evaluated using the 
CTEQ6M PDF set with all scales equal to 50 GeV. The errors represent Monte Carlo statistics 
only. Note that the values of Ft at LO and NLO are 1.651 GeV and 1.480 GeV respectively and 
the branching ratio of the W into leptons is Br (IF — ► ev) = 0.1105. 



Collider, y/s 


aoB t -,bue [fb] 


oBt^bve [fb] 


aBt^bue+x [fb] 


pp, 1.96 TeV 


0.8564 ± 0.0006 


0.7887 ± 0.0005 


0.7806 ± 0.0005 


pp, 14 TeV 


356.9 ± 0.2 


391.7 ± 0.3 


395.7 ± 0.3 



The results shown in the table agree with this expectation. Since, in our approach, the effect 
of the NLO corrections in the production is fairly small, the numerical difference is only 1% 
at both colliders. 

We conclude this section with a more detailed presentation of the scale dependence of our 
calculation when using the updated parameters and PDF set. Anticipating the study of the 
following section, we also choose a lower value for the b jet veto in the next-to-leading order 
calculation, veto = 30 GeV. This results in the scale dependence shown in Fig. where 
the cross sections at each choice of scale are expressed as a ratio with the central result at 
fi = veto and we vary the scales by a factor of two about Pq. We also show the curves 
obtained when varying the renormalization and factorization scales separately. One sees 
that the relatively small scale dependence at lowest order is the result of a large cancellation 
between the dependence on the factorization and renormalization scales individually. In 
contrast, the dependence on these scales at next-to- leading order, either on their own or 
when varied together, is small - less than 10% over this range. We also note that, for this 
choice of parameters and veto, the next-to-leading order corrections do not alter the tree 
level cross section, which remains at 346 fb. 
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FIG. 9: Scale dependence of the cross sections for W~{— ► e~v)t{— > ue + b) production at the LHC, at 
LO (left) and NLO (right). The NLO calculation includes the effect of radiation in the decay of the 
top quark. The scale /i is expressed as the ratio with the central scale /io = Vt veto = 30 GeV and 
the cross sections are scaled to the central result at \x = /io- The solid lines represent the variation 
of both renormalization and factorization scales together (fiR = [if = n), the dashed ones the 
result when only hf is varied (fiR = /io) and the dot-dashed curves represent the dependence on 
the renormalization scale alone {^F = Ho)- 



V. THE Wt BACKGROUND TO H -> WW* 



As an example of the utility of our calculation, in this section we consider the effect 
of NLO corrections in the context of a search for the Higgs boson at the LHC. In the 
intermediate mass range, 155 < mjj < 180 GeV, one of the search strategies involves Higgs 
production via gluon fusion, with the subsequent decay of the Higgs boson into off-shell W 
pairs which then decay leptonically 



g + g^H -> W~ +W+ 



(14) 



e + v 



The largest background in this channel is from the continuum production of W pairs, both 



from di 
fusion 



joson production via quark- ant iquark scattering and from loop-induced gluon-gluon 
28l l29j | . A further significant source of background events comes from processes 



producing top-quarks that decay leptonically. Since the presence of neutrinos in the signal 
prevents a full reconstruction of the Higgs mass peak, an accurate prediction of all the 
backgrounds is necessary. 
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In previous studies [30, 31] two top backgrounds have been considered using Pythia [32]. 
These are resonant it production and Wt production, the process that we consider here. 
As we have discussed previously, these two processes become entangled at NLO. We will 
separate them according to the procedure that we outlined in Section 11111 Therefore, in 
addition to our NLO calculation of the Wt process, we will also consider the contribution 
from resonant ti production, with the appropriate top quark decays. 

Since the signal process contains no jets at lowest order, it is efficient to impose a veto 
on all jet activity to reduce the size of these backgrounds. Our application of a cut on the 
transverse momentum of the b jet in the Wt process fits naturally into this procedure. We 
simply extend our veto to disallow all contributions with any jet observed above the veto 
threshold. We note that there is a slight mismatch due to the fact that our theoretically- 
motivated veto applies at all rapidity values whilst the experimental approach only vetoes a 
jet up to a few units of rapidity. However, we do not expect this to greatly affect our results. 

For our parton-level study we adopt a minimal set of cuts and examine the effect of 
the NLO corrections on a selection of observables that are typically used in more detailed 
expe— studies HHHQ. Ur basic cuts represent the unite acceptance ■nn.ts 
of the detectors at the LHC, 

p r (lepton) > 20 GeV, |r/(lepton)| < 2.5, (15) 

applied to both of the leptons produced in the W decays, with the missing transverse 
momentum also constrained by, 

Pr (missing) > 30 GeV. (16) 

The final cut that we apply is the jet veto, after potential jets have been clustered according 
to the &t algorithm with a jet separation parameter AR = 1.0. Events are not included if 
any jet is observed with, 

pr(jet) > 30 GeV, |r/(jet)| < 3. (17) 

We note that, in addition to excluding additional radiation at next-to-leading order, this 
veto also applies to the b jet that is produced in the top quark decay. This results in a 
substantial decrease in cross section compared to the totally inclusive case. 



To exploit the spin correlation between the leptons in the signal events 



33j , one can make 



quite stringent cuts on the opening angle between the leptons in the transverse plane, A(/)£ 
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FIG. 10: The distribution of the opening angle between the leptons in the transverse plane, 
for the signal and the two background processes considered here. All curves are lowest order 
predictions and are normalized to unity. The signal calculation uses a Higgs mass of 155 GeV and 
the backgrounds are W~t (solid) and tt (dashed). 

This is illustrated in Fig. where we show the shapes of the lowest order predictions for 
the Higgs signal and the two top backgrounds, Wt and ti. Signal events predominantly 
contain leptons with a small opening angle between them, whereas both backgrounds tend 
to produce leptons that are almost back-to-back in the transverse plane. 

We now examine the extent to which this is changed in our NLO calculation of the Wt 
background, with our results shown in Fig. ^2 The effect of the NLO corrections is to 
change the shape of the distribution considerably when these cuts are applied. The peak at 
large is shifted to a smaller value and becomes much less pronounced. This could have 
quite a large impact on a strategy in which this background is measured using events at 
large A0« and then extrapolated via the theoretical shape to the Higgs signal region. One 
also sees that the shape is changed again when including the effects of radiation in the top 
quark decay, with the peak being sharpened once more, although the effect is fairly minor. 

One can also imagine constructing the transverse mass of the putative Higgs boson (mr) 
from the transverse momenta of the dilepton system and the missing p?, 




(18) 
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FIG. 11: The distribution of the opening angle between the leptons in the transverse plane for 
W~t events at the LHC, calculated at LO (dashed) and NLO (solid) for the LHC. The NLO curves 
are labelled according to whether or not they include the effect of radiation in the decay. All the 
rates are normalized to unity. 

where A</> is the angle between the two vectors in the transverse plane. Cutting in a suitable 
mass window can further help to reduce the backgrounds for only a small loss in signal. The 
impact of our next-to-leading order calculation on this distribution is illustrated in Fig. 
In this figure we have only shown the result when including radiation in the decay but note 
that this distribution changes little when it is excluded. One can see that the shape of this 
distribution is relatively unchanged at NLO, although more events are produced at high 
values of M T , beyond the peak of the distribution, than at LO. 

Finally, to give some idea of the effect of the NLO corrections on the number of events 
that should be observed in this channel, in Table ILTT1 we show the cross sections that we find 
for the Wt and ti processes. The W's decay into electrons only and for the Wt process, both 
W~t and W + i are included. Results are shown at LO, NLO and at NLO when including 
radiation in the top quark decay. For the ti process, the lowest order diagrams of Fig. |H1 (b) 
are calculated. As an approximation to the NLO result, a i^-factor is applied from the NLO 
calculation involving no top quark decay. In order to match the study more closely, we have 
used a i^-factor obtained when applying the jet veto of Eq. (|17|). Using the common scales 
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FIG. 12: The distribution of the transverse mass for W~t events at the LHC, calculated at LO 
(dashed) and NLO (solid) for the LHC. The NLO result includes the effect of radiation in the 
decay. The rates are normalized to unity. 

fjjji — fMp — rrit, we find that this factor is /C = 0.7. 

The cross sections shown in Table ITTT1 are evaluated using three different sets of cuts. The 
first column uses just the basic set of cuts (Eqs. (|15j) - (fT7|) ). then the other two columns 
represent extensions of these cuts that might be used in the search for a Higgs boson of 
given mass. Both the further sets also impose, 



to select the Higgs signal region. In addition we have used a cut on the reconstructed 
transverse mass around the Jacobian peak of the putative Higgs mass. In the first case, the 
cut is constructed for a Higgs mass of 155 GeV, by constraining 125 < < 155 GeV. The 
second set requires that 140 < tut < 180 GeV and is aimed at a search for a 180 GeV Higgs. 

One sees that the effect of the NLO corrections in this region of phase space is significantly 
different from the inclusive case. The cross section is decreased substantially, with a i^-factor 
(for the calculation including radiation in the decay) of approximately 0.6 when applying 
the Higgs search cuts. We also see that the contributions from the two top quark processes 



A(p u < tt/4, 



vau < 35 GeV, 



(19) 
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TABLE III: Comparison of LO and NLO cross sections for top quark backgrounds in an inter- 
mediate mass Higgs search at the LHC. Results are shown in femtobarns for (W~t + and 
tt production, with leptonic decays of both the W and the top quark. Three sets of cuts are 
considered, which are described in detail in the text. The NLO Wt calculation is performed both 
without including QCD effects in the decay (aB t ~,bve) and also when it is included (aB t ^bve+x)- 





basic 


m H 


= 155 GeV 


m H = 180 GeV 






Wt process 






<7 B t ^bue [fb] 


40.08 




0.80 


0.80 


crB t ^ bue [fb] 


13.14 




0.38 


0.43 


cBt^bve+X [fb] 


13.85 




0.46 


0.51 


a LO x K. [fb] 


30.52 


tt process 


0.42 


0.43 



are comparable when calculated in this way. 

VI. CONCLUSIONS 

At the LHC, the top quark will be readily produced in association with a W boson. In 
this paper we have performed a next-to-leading order calculation of this process, including 
both the subsequent leptonic decays W~ — > e~v and t — > ue + b, as well as the emission of 
real radiation in the top decay. 

For total inclusive cross sections, where the top quark and W~ boson do not decay, 
comparison with previous calculations is possible and we find results which are broadly 
similar. However, due to the presence of doubly-resonant tt diagrams at NLO - which 
are handled differently in our calculation - we do not find exact agreement. Our method 
maintains the consistency of the 6-PDF approach and requires the simultaneous use of a 
relatively low factorization scale hf ~ ijn t + mw) / ^ and a veto on b quarks with a transverse 
momentum larger than this value of hf- 

With this approach, we find that the NLO corrections to the inclusive rate at the Tevatron 
can be large, but the process remains phenomenologically irrelevant there. At the LHC 
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the corrections are smaller and can either increase or decrease the cross section by 10%, 
depending on the b quark veto that is chosen. We also performed an analysis of this process 
as a background to Higgs production at the LHC, where quite a severe veto on all jet activity 
is applied. We find that the shape of the A0« distribution - crucial to the estimation of 
this background - is significantly changed. However the cross section for this process can be 
reduced by as much as 40% in the Higgs signal region. 

Our predictions are based on a number of approximations. Firstly, the mass of the b quark 
is set to zero throughout and the top quark is kept on its mass shell. The inclusion of the 
6-mass in the top quark decay and implementation of a Breit-Wigner distribution for the top 
quark can be studied in the lowest order calculation. Neither of these significantly alters the 
total cross section or the shapes of the distributions which we have examined. Effects due to 
the interference between radiation in the production and decay stages are estimated to be 
similarly small, of order a s T t /m t . Lastly, our predictions are performed at parton-level only 
and lack any modelling of hadronization and showering. This is particularly relevant to our 
phenomenological analysis, where the application of a jet veto at relatively low transverse 
momenta (compared to the top quark mass) increases the sensivity to soft-gluon effects. 
Nevertheless, we hope that our results can provide a starting point for further investigation 
of the Wt process at next-to-leading order accuracy. 
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APPENDIX A: INTEGRATION OF DIPOLES 

We have slightly extended the approach outlined in Appendix A of Ref. y in order to 
handle the presence of gluons in the initial state. We refer the reader to Ref. 101 and the 
original papers CS [13] and CDST |14| for further explanation of the method and notation. 
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1. Initial-state emitter with initial state spectator 

As explained in Ref. [l(|, we have generalized the dipole phase space (CS, Eq. (5.151)) 
by introducing a constraint enforced by the factor 0(a — The variable Vi is the rescaled 
value of the propagator defined by, 

PaPi / a -i \ 

Vi = (Al) 

PaPb 

where p a is the initial state emitter, p^ is the emitted parton and p b is the other initial 
state parton which is the spectator. Further details are given in section 5.5 of CS. In this 
appendix we extend the previous treatment by considering the remaining q, g and g, g cases. 
The dipole integrands which we subtract are obtained by modifying CS, Eq. (5.154): 



(V M (x hab )) = 8n^a s C F 



l + a-Xi — 



'H,ab) 



•Ei,ab 



~V € x i,ab 



(A2) 



_ 1 •Kijdb %i,ab 

We find that the result for the q, g case is given by, 

l + (l-x) 

X 



•Ei,ab 1 %i,ab . \ 

T r Xi,ab \± •^i,abj 



V q ' g {x;e,a) = C F 

1 fa 

2 ln(l - x) - - + 6(1 - x - a) In 

e \ 1 — x 



+ r]x } + 0(e) 



and for the g, g case it is, 



V 9 > 9 (x:e,a) = C, 



7 ~ ~6 



6(1 -x) 



+ 2(x(l-x) + - — - - l) ^2 ln(l - x) - - 

. 2(x(l -x) - l) 2 , / a 
+ 9(1 - x - a) v v - — '— — -In 



[1-xl 



+ 4 



x(l — x) 

ln(l — x) 

1 — x 



.1 — x 
+ 0(e). 



(A3) 



(A4) 



(A5) 



In the limit a = 1 these functions correspond to those given in CS, Eq. (5.155) and also 

nn 

agree with the results of Nagy [34, 35]. 
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2. Initial-state emitter with final-state spectator 



In this category, we complete the treatment of Ref. [10] by considering the remaining 
three cases: 

a) Initial g — > q + q, 

b) Initial q — > g + q, 

c) Initial g — ► g + g, (A6) 

where the first parton on the right-hand side is an initial-state emitter and the last one is 
a massive final-state spectator. The phase space is the generalization of CDST, Eq. (5.79) 
with an extra factor of Q(a — Zi). 

For case (a) the dipole integrand is given by a generalization of CDST, Eq. (5.82). The 
result is written in the terms of the variable z + defined by, 

1 — x 



2+ 



2 • 



1 - x + n Q 



(A7) 



and we find, 



I 9q {x- e,a) = T R Ua? + (l-x) 



x 



In 



:i-x) 



1 — X + 2/1 2 



--6(z + -a)ln ( — 
e \ a 



In this expression we have also introduced the variable, 



~ 2 V 



m 



+ 2r]x(l-x) } +0(e) . (A8) 



(A9) 



x 2p ai p j 

which only depends on p ai and pj (defined in CDST Eq. (5.73)), the momenta held fixed 
when the x integration is performed. 

For case (b) the integrand is obtained by generalizing Eq. (5.84) of CDST. Performing 
the integrals yields the result, 

Iq [x] e,a) = C F < 



x 



X 



111 



(1-x) 2 
1 — x + xfl 2 



Q(z + — a) In ' 



a 



2p 2 In 



~ 2 
Xfl 



1 — x + xfl 



r,2 



Q{z+ - a)2fi 2 In ± 

V 1 — a 



+ r]x } +0(e) 



(A10) 



24 



Lastly, the dipole integrand for case (c) is an extension of CDST, Eq. (5.86). The result 
in this case is, 



I 9 Q 9 (x;e,a) = C A l5{l-x) 



- + - ln(l + /i 2 ) + — 



e 2 e 



1 



+ 2 Li 2 (-/i 2 ) + 2 ln(/} 2 ) ln(l + /i 2 ) - ± ln 2 (l + /i 2 ) 



2 Mi+^) +4 



e [1 — x]+ [1 — x 

1-x 



ln(l — x) 



1-x 



+ 



+ 2 (x(l -x) + - l) - i + In ( J^L) - 6„ + - „) .„ (it) 



— @(-2+ — a) 



2jl 2 In - In -( 

V 1 - a / (1 — x) — x + a)/ 



+ 2 ^l n (-^)--^l„( 
VI — x + xr/ 1 — x V 



X + XfJ, 



2 

(1-x) 



(1-x) Vz+( 
2 — x + x/i 2 



1 + /2 2 



)}■ 



(All) 



APPENDIX B: THE VIRTUAL ++ AMPLITUDE 



We remind the reader that we calculate the virtual amplitudes for the production of a 
W that decays leptonically in association with an on-shell top quark, 



b + g — > W(-^ e + n)+t, 

where n represents the neutrino, t 2 = m 2 and b 2 = g 2 = e 2 = n 2 = 0. 
We first introduce the following short-hand notation, 



(Bl) 



<xy 



= 2 x ■ y , txy = (x + y) 2 , Q = q 2 = 2 



e • n 



Qxy — Q txy j Q — Q Tn>t 

S tb = (bt ) [t b] = 2 b ■ t = 2 b ■ t - m\ 



b ■ g 



(B2) 



t 4. ' 

t-g 

such that the massive top quark momentum t does not appear directly in the following 
formulae but is instead replaced by the massless vector to, which is suitable for our spinor 
approach. 

Next we define a set of functions that are useful for decomposing the form of the virtual 
amplitudes. In these formulae, the scalar n— point functions are written using the following 
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definition, 

TOPi...n-i}; {<.„}) = ^ IS 

r T J (2tt 

where, 



)4-2e {l 2 _ _ pi) 2 _ ^ ((/ _ ^2 _ <} > 

(B3) 



r r = 



mr 



2\ e 



16 7r 2 r(l -e)' 



We subsequently identify J 2 (. . .) = B (. . .), / 3 (- • •) = C {. . .) and J 4 (. . .) = D 
simplest functions contain only bubble integrals, 



Q = l/e + 2-B (q;0,m 2 )=Qlog(-Q/m 2 )/Q 



(B4) 



The 



= 1/e + 2 - £ (* + 0; 0, m 2 ) = r tg \og(-r tg /m 2 )/t 



tg 



= 1/e + 2 - B (t + b- 0, m 2 ) = r tb \og(-r tb /m 2 )/t tb 

= Q (kg ~ Iq) I 'qtg 

= t tb (l tb — lo)/qtb 

= (kg - h) Q 2 /q 2 tg + h Q 2 /qtg/n g - Iq Q/qtg + Q/qtg 



(B5) 



Three triangle functions appear in our results, which we choose to keep as independent 
functions, 



= C (t + b,g;0,m 2 ,m 2 ) 
= C (t,g;0,m 2 ,m 2 ) 
= C (g + b,t;0,0,m 2 ) 



(B6) 



Lastly, we form functions which are combinations of the triangle functions and the basic box 
integrals that enter our calculation, 

L\ = T tb T gb (D (g,b,t;0,0,0,m 2 )-C (g,b;0,0,0)/r tb 

+C (g,t + b;0,0,m 2 ) (r gb + r tg ) / r gb / r tb - C (b, f; 0, 0, m 2 )/r gb 

+C (g + b, t; 0, 0, m 2 ) (r tb (r tg + r tb ) - 2 m 2 r gb )/r gb /r 2 b ), (B7) 



r 6 

L 2 



TtbTtg (D (b,t,g;0,0,m 2 ,m 2 ) - C (b,t;0,0,m 2 )/T tg 

+C (b, t + g;0, 0, m 2 ) (r gb + T tb )/r tg /T tb + C (t, g; 0, m 2 , m 2 ) (2 m 2 r gb - r tg T tb )/T tg /T t 



+C (t + b,g; 0,m 2 ,m 2 ) (r tg + r gb ) (r tg T tb -2m 2 T gb )/T 2 g /r 2 b ), 
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(B8) 



L 6 3 = T tg T gb (D (b,g,t;0,0,0,m 2 t )-C (b,g;0,0,0)/r tg 

+C (b, t + g; 0,0, m\) (r gb + T tb )/r gb / r tg - C (g, t; 0, 0, m 2 t )/T gb 

+C (g + b, t; 0, 0, m 2 t ) (r tg (r tg + r tb ) - 2 m\ r gb )/r gb /r 2 g ). (B9) 



All the basic scalar integrals, Bq, Cq and Dq are well known [36, 37]. 

To renormalize the virtual amplitudes we have used the modification of the MS scheme 
in which the top quark is decoupled. The top self-energy is renormalized on-shell so that 
we have evaluated the self energy in the top- right diagram of Fig. |2] including the following 
mass counterterm, 



«W = -£ (±f) j^y (t + 5-i| **, (Bio) 



e 



and we renormalize the top wave function by adding, 
SZ w f a s /4 7r/i 2 V Cp ( 3 



2 4tt V mi J 2T 1 - e 



e 



- + 5 - T]) A Tree Leve i. (Bll) 



In the above formulae rj is a parameter that specifies the regularization scheme adopted 
throughout the calculation. 77 = corresponds to the 4-dimensional helicity scheme, which 
we use here, whilst rj = 1 is appropriate in the 't Hooft Veltman scheme. Subtracting 
the top loop contribution to the gluon self energy at zero momentum transfer enables the 
coupling constant to evolve due to the presence of 5 light flavours only. Finally, the partial 
cancellation of the coupling constant renormalization and this top contribution to the gluon 
self energy give us the last contribution, 



SZ ' g = ^ (4 ,). " W / 6 -y A Tm M . (B12) 
47T el 1 — e 



We are now in a position to write down the amplitudes for the specific helicity choice 
'++'. The first '+' signifies the helicity of the gluon and the second '+' the spin of the top 
quark in a basis determined by our decomposition in Eq. (jB2|) . The tree-level amplitude is, 

A ++ _ V^9s9wT a ++ 

A Tree Level ~ Ty_ ^0 > \ Dlc> ) 

prop 



where, 



and, 



M + = ~m t ^^"^y «bto) [bn] + (gt ) [gn]), (B14) 



W- rop = q 2 -M 2 v + iM w Y w . (B15) 
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The corresponding virtual amplitude is decomposed into a piece containing poles in e which 
is proportional to the lowest order result, plus a finite remainder, 

ku^ - ^i^S^f)' [(F+p ™> yAt++At+] - (B16) 

The pole pieces are given by, 

F = -N (3/2/e 2 - 1/e log(-r tfl / m 2 ) - 1/e \og(-r gb /m 2 t ) 
+l/2/e + log(-r tg /m 2 ) 2 + 1/2 \og(-T gb /m 2 ) 2 ) 
+1/2/ N (1/e 2 - 2/e log(-r t6 /m 2 ) + 1/e + 2 log(-r t6 /m 2 ) 2 ) 
-N (U 2 (t tg /m 2 ) + tt/12) + (Li 2 (^/m 2 ) + 7r 2 /12)/iV, (B17) 

while the term which is the sum of coupling constant and wave-function renormalization is, 

P wf = (2T R n f /3 - 11 N/6) (1/e - log^/m 2 )) - 3C F /2 (1/e + 5/3). (B18) 

The remainder is then written in terms of our functions as, 

At+ = wwwU] & [L ' :s,i/2 + L% {Sa/2 - Mb - q " m ' /VT " ] 

-Cq r gb 5 tb m 2 / r tg -2Cq 5 tb m 2 + 2 Cjf 5 tb m 2 + C$ r gb 5 tb m 2 / r tb + l] b r gb t tb / r tg + l Q 5 tb 

—kg Ug Stb/ T tg — kb ttb Ttb/ T t g + kb T gb Stb/ Ttb + Ttb W T tg] + Cp [l tg q tb t tg / T tg — l tg Q — 5 tb 

+ l t g Sfb + 2 l tg t t g 5 tb / Tfg + Qtg fa/ Ttg ~ 2 S tb Q / Tfg] +N[- LI 5 tb / 2~Cq Tg b 5 tb 7^/ Ttg] } 

+ W(lt) 2 \bto] i ^ ^ Ttb/2 + L " ^ ~ ttb/2 + ^ Ttg,Tgb) + 2 °° ^ Tgb 5t " m2t,Tth/Tt9 
+Cq q tb S tb m 2 / r tb - Cq q tb S tb m 2 / r tb - Cq r gb S tb m 2 / r tb -2C$ r gb 5 tb m 2 / r tg - l tb T tg S tb / r tb 

~k g S t b\ + C F [ - 4 Cq r gb Stb m 2 / T tg + 21q r gb Q/ r tg + l Q r gb - 2 l tg q tb - 2 kg qtb Ql n g 

+k g T tg Q/t tg + Iq q tg qtb/ T tg + l Q Q qtb/ T tg + Iq q t b + IqQ~ Qtg r tg / t tg + r gb - S tb ] 

+N [Ltr tb /2-LtS tb ]} 

+ \gbHgt )\bC)\bt } { ~N~ I ~ L2 qtb m ^ 2 / Ttg + °® q ' b 6tb m * ^ Ttg ~ °° qtb 6tb m ^ Ttb 

-Iq T gb Ql T tg + kb Q Stb/ T tb -IqQ qtb/ T tg + l t g q t b Ql Ttg] + C F [l 2 g TgJ T tg + 2 Iq T 2 b Q/ T? g 

-3 l Q r gb/ T tg - 2l Q T gb + 4 kg Qtb Q/ T tg ~ 2 l tg q 2 b Ql T 2 g - 2 kg t tg q 2 J T 2 g + 2 kg q 2 J T tg 
-kg T tg Ql t t g + Iq q tg ql m 2 / t 2 J Q + l Q Q q tg q\ b / t 2 J Q - 2 l Q Q q tb / r tg + Iq q 2 b m 2 / T tg / Q 
+h Q Qtb/ Ql r t g ~ qtb Stb/ T tg - q 2 J T tg + T tg Q /Ug] Y (B19) 
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The expressions for the other amplitudes are similar but slightly longer, requiring the addi- 
tion of a further 5 functions to describe them compactly. We do not reproduce them here, 
but they are available as Fortran files from the authors on request. 
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At the Tevatron, this is due to the fact that the machine is a proton anti-proton collider. 
At the LHC, the equality is because the perturbatively-derived bottom quark distribution 
functions that we use are the same for b and b quarks. 

This can be seen by examining a study of the similar process, charged Higgs production in 
association with a top quark [19j. We have also reproduced this result for the process at hand. 
A NLO calculation including the decays of the top quarks is not currently available. However 
one could, for instance, use a lowest order calculation including the decays and normalize to 
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the inclusive tt NLO rate, which is known 



2C|. 



As we have already pointed out, interference effects between singly- and doubly-resonant 
diagrams are small. Although this procedure is not strictly gauge invariant, it is no more serious 
an error than that incurred when introducing a Breit-Wigner width for resonant propagators. 
We use CTEQ5L1 for the lowest order calculation and CTEQ5M1 at NLO, the versions that 
include the improved evolution code. 
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